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We consider the equilibrium statistical properties of interfaces submitted to competing interac- 
tions; a long-range repulsive Coulomb interaction inherent to the charged interface and a short-range, 
anisotropic, attractive one due to either elasticity or confinement. We focus on one-dimensional in- 
terfaces such as strings. Model systems considered for applications are mainly aggregates of solitons 
in polyacetylene and other charge density wave systems, domain lines in uniaxial ferroelectrics and 
the stripe phase of oxides. At zero temperature, we find a shape instability which lead, via phase 
transitions, to tilted phases. Depending on the regime, elastic or confinement, the order of the 
zero-temperature transition changes. Thermal fluctuations lead to a pure Coulomb roughening of 
the string, in addition to the usual one, and to the presence of angular kinks. We suggest that such 
instabilities might explain the tilting of stripes in cuprate oxides. The 3D problem of the charged 
wall is also analyzed. The latter experiences instabilities towards various tilted phases separated by 
a tricritical point in the elastic regime. In the confinement regime, the increase of dimensionality 
favors either the melting of the wall into a Wigner crystal of its constituent charges or a strongly 
inclined wall which might have been observed in nickelate oxides. 
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I. INTRODUCTION 

Various types of systems display peculiar properties 
due to competing long-range forces. In this paper we fo- 
cus on the statistical properties of charged interfaces, in 
Ising-like systems, such as aggregates of charged topolog- 
ical defects in charge-density wave systems^** and among 
them polyacetylene^, charged domain-lines in uniaxial 
ferroelectrics* or stripes in oxides 5,6 . Instabilities, such 
as the observed inclination of stripes in cuprates* and 
manganese 8 , might be related to the presence of the long 
range Coulomb interactions and its competition with an 
attractive force. Uniaxial ferroelectrics and density waves 
are also model systems where such instabilities could 
take place. The present study deals with one or two- 
dimensional charged interfaces, i.e. strings or walls. The 
former is realized in oxides, monolayers of doped conduct- 
ing polymers or other charge-density waves and junctions 
in field effect experiments* for equivalent materials. The 
latter might be present, due to a dimensional crossover 
which induces a 3-D ordering, in the systems cited above. 
The above interests emerged from a model-independent 
theory concerning the statistical and thermodynamic 
properties of uncharged and charged solitons which has 
been studied, respectively, in Refs. ITol and fill The sys- 
tem considered is two-fold degenerate and, when topolog- 
ically doped, the solitons form a one-component plasma 
in either 2 or 3D space. The competition between the 
long-range 3D Coulomb interaction of the particles and 
the confinement force between them has led to a very rich 
phase diagram. Particularly interesting was the regime 
where the temperature is much less than the confine- 
ment energy scale, so that the solitons are actually bound 
into pairs, and the repulsive Coulomb interaction is weak 
enough to preserve these bisolitons as the elementary par- 
ticles. In this regime, aggregated phases of bisolitons 
were shown to exist. The study of such aggregates in 
the continuum limit is equivalent to the study of charged 



interfaces. 

The statistical properties of charged fluctuating mani- 
folds have been considered before!*. In this study, the 
authors have been considering a generalization of the 
theory of polyelectrolytes, cf. Ref. 1 1.1 and have been fo- 
cusing on the properties emerging from scale invariance 
and renormalization. Even though some of our results 
will converge with the previous study, which seems to 
be quite natural, our aim will be different. We consider 
interfaces directed along an anisotropy axis of the sys- 
tem which are either in an elastic or confinement regime; 
the latter is relevant to quasi-one dimensional systems in 
relation with the statistical properties of solitons men- 
tioned above. In both regimes the Coulomb interaction 
favors the disintegration of the interface. In contrast with 
the instability towards a modulation which has been sug- 
gested in Ref. ITU this competition, as we will show, leads 
to zero-temperature shape instabilities with respect to a 
tilting of the interface. We will focus mainly on the prop- 
erties related to this discrete symmetry breaking. This 
will able us to better understand the structure of the 
solitonic lattice in the presence of the long range repul- 
sive interaction. In this respect, and as explained inii, 
even though this problem seems to be quite old, the de- 
tailed impact of the Coulomb interaction has not been 
reported before. Moreover, the charged solitonic lattice 
is commonly observed nowadays as stripes in various ox- 
ides. Experimentally the observation that upon doping 
the system, at low temperatures, a transition from a 
collinear, i.e. vertical, to a diagonal stripe phase takes 
place* is strikingly similar to our actual results. 

The paper is organized as follows. We will mainly con- 
centrate on the statistical properties of a charged string. 
In section ITT1 we present the models of the confined and 
elastic charged strings. In section ITTTI we derive the sta- 
tistical properties of the string with the help of a saddle- 
point calculation. We present the zero-temperature re- 
sults dealing with the instabilities of the string. In sec- 
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tiou lIVI the effects of thermal fluctuations are considered. 
In section IVI the previous results are confirmed with the 
help of a numerical approach. In section IVTI the present 
study is applied to the 2D and 3D solitonic lattices, as 
well as stripes and uniaxial ferroelectrics. 

II. ELASTIC AND CONFINEMENT REGIMES 
FOR THE CHARGED STRING 

The problem of the charged stringil has emerged from 
the study of the statistical properties of topological de- 
fects, i.e. solitons, common to charge density wave 
(CDW) systems. These general arguments are repro- 
duced in the following as they give a full meaning to the 
physical origins of the present theory, cfi for a review 
on the physics of topological defects in CDW systems. It 
should however be clear that this study is also relevant 
to other systems with discrete symmetry breaking as has 
been mentioned already in the Introduction. 

For the moment we consider a one-dimensional system, 
i.e. a single chain, in a CDW state. The latter is de- 
scribed by a lattice order parameter which is complex in 
general 

A (a;) = |A(a;)|exp(i<p(a;)) 

where x runs along the chain, |A(x)| is related to the 
CDW gap and (f(x) fixes the position of the CDW with 
respect to a host lattice. In weakly commensurate CDW 
the low energy excitations are excitations of the phase. 
This leads, cf^, to a sine-Gordon type Hamiltonian 

*-/*(£(«)'-" r ""H (i) 

where the first term gives the elastic energy of the CDW 
and the second term reflects the pinning of the CDW 
superstructure by the host lattice; M is the degeneracy 
of the ground-state. 

The non-trivial solutions of the equation of motion asso- 
ciated with (JJJ are the solitons, i.e. the (p— particlesi*, 
which may be viewed as compressions or dilatations of 
the CDW. As we are in the realm of electronic crystals 
the solitons might therefore carry a charge which is frac- 
tional, i.e. q = 2/M, in the general case. For a double- 
degenerate ground state, M = 2, they are given by 

ip(x) = 2 arctan (exp(2x/£)) (2) 

where the length £ ~ y/Hvp/W corresponds to the width 
of this nucleus domain wall connecting the two ground 
states ip(±oo) — ±1. 

The situation is slightly different for the two-fold degen- 
erate trans-polyacetylene because in this case the band is 
half-filled and the system is commensurate. In particular 
equation |QJ does not describe this case. However the 
previous arguments remain qualitatively the same. The 



defect is an amplitude or n— soliton; it carries an integer 
charge and is described by 

A(x) = A tanh(a;/0 (3) 

where Ao is a constant deformation of the chains. Again 
this excitation connects the two ground states of the sys- 
tem A(±oo) = ±A . 

For a quasi-one-dimensional system, which we consider 
here, neighboring chains must have the same phase mod- 
ulo 2ir. This leads to the following general two-particle 
interaction 

E(5x) = -V J dzA a (Sx - z)Ap(z) (4) 

between neighboring chains a and (3, with coupling con- 
stant V , and solitons distant by Sx. In the coarse-grained 
regime, Sx > £, © grows linearly, E » VSx, a signature 
of the confinement of the defects. This can be shown 
explicitly with the help of © and @. Dealing with a 
string, e.g. a domain line of solitons, this confinement 
energy is equivalently reproduced by the Solid on Solid 
(SOS) model. The related Hamiltonian reads 

j L-l 

Hconf = — - X v\ ( 5 ) 

» J/=0 

where we have taken J± = Va y as the inter-chain energy 
scale, y runs along the vertical direction and x y is an 
integer variable giving the deviation of the string with 
respect to the y axis. 

The opposite case, corresponding to small deviations 
Sx <C £, is more general for the energy then increases 
quadratically E w VSx 2 . This is equivalent to the elastic 
model which is extensively used in interface physics. The 
related Hamiltonian reads 

H elast = J ± [^(^) 2 . (6) 
J a y \dyj 

These two interface models are well known 16 . To com- 
pactify notations we shall define 

J % ay 

where p = 1,2 leading respectively to the confinement 
and elastic regimes and where the continuum limit has 
been taken in both cases. We thus see that a crossover 
(corresponding to fractional 1 < p < 2) from the elastic 
to the confinement regime takes place as the distance 
between solitons increases, i.e. as the tilt angle of the 
string increases. This is manifested by the fact that Sx = 
(dx/dy)a y = Qa y where ( is related to the tilt angle of the 
string and a y is the inter-chain distance. For £ <C £,/a y 
the string is elastic whereas for ( ^> £/a y the confined 
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string is relevant. The crossover between the two regimes 
takes place at 

CO = i/Oy. (8) 



where z < 1 so that the elementary constituent of the 
string might carry a fractional charge ze (z = 1 for poly- 
acetylene), e is the dielectric constant of the isotropic, 
neutral media in which the plane is embedded and the 



III. THE CHARGED STRING AT ZERO 
TEMPERATURE 

A. The saddle point approximation 

We consider a configuration 

x(y) = x (0 \y) + 6x(y) (11) 



where H p is given by GDJ and H p 0) = H p {Sx^(y)}. The 
saddle point distribution is given by SH p /Sx(y) = that 



This non-linear integro-differential equation cannot be 
solved exactly. We will therefore take an ansatz for the 
ground-state configuration a;^ ^ (y) of the string. As can 



The introduction of the 31? long-range Coulomb interac- 
tion yields 



(9) 

I 

contribution of the vertical string has been subtracted. 

The thermodynamics of the charged interface is governed 
by H p = H 0p + H c or explicitly 



(10) 



were x{y) has been expanded in the vicinity of the sad- 
dle point distribution a;' ' [y] with deviations 6x(y). Free 
boundary conditions are taken. The expanded Hamilto- 
nian reads, in the harmonic approximation 



(12) 



is 



(13) 



I 

be seen from equation l|10|l . at zero temperature, the elas- 
tic or confinement term favors a vertical string, i.e. along 
the y axis, whereas the Coulomb favors the disintegration 
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of the string in the perpendicular direction. Instabilities 
of the string arise from the competition between these in- 
teractions. We find that these instabilities lead to a tilted 
string as the new ground state of the system. Therefore 



(14) 



where £ = tan#, 9 being the tilt angle of the string with 
respect to the y axis. This ansatz satisfies the linearized 
equation l|13l) up to some logarithmic corrections. As will 
be shown explicitly in the following, £ is related to the 
ratio between the Coulomb and confinement or elastic 
energy scales. It is thus related to the elementary charge 
of the string constituent soliton. 



B. The elastic charged string 

We first consider the elastic regime, with p = 2 in 119L 
corresponding to tilt angles below the crossover value JSJ. 
Figure ^ displays this free energy as a function of the 
angle. Until 



7d = 2 



(20) 



the vertical line, £ 2 = 0, is stable. Above 7^ the vertical 
line becomes unstable and the free energy has a double 
well shape. This is due to the double degeneracy of the 
system; the new, tilted, ground states corresponding to 
±£2 have the same energies. The optimal angle reads 



With the help of JHJ, and (HHJ, (EJ reads 

k 

with the eigenvalues 



(15) 



4AC P ) = k 2 Ji_ \p(p-l) + 



Kl - 1 



[1 + C P 2 ] 5 



r7(*0 • (16) 



Q p is the optimal angle corresponding to the distribution 
minimizing H p 0) and 



with 



7(fe) = 7log(l/fca) 



w y (ze) 2 
1 = JI U * = 2^' 



(17) 



(18) 



7 being the ratio of the Coulomb energy scale w y to the 
elastic or confinement energy scales. 
The first term in i|15fl corresponds to the mean-field 
Hamiltonian. Hp ' 1 leads to the following mean-field free 
energy density 



/f (0 = - 

a y 



C P + 7D 



- 1 



(19) 



^2 



(?)' 



(21) 



The transition at j D is of the second order. This can 
be shown by evaluating a quantity analogous to the heat 
capacitance in usual thermodynamic phase transitions. 
That is 



C 7 = -JD 



(22) 



with the tilt angle given by l|21|l above the transition and 
vanishing below. Eq. 122|) yields 



2Jl 
~9~ 



Id 



2/3 



i(y)" ! 



7d > Id 



and C 7 = otherwise, which leads to a jump at the 
transition, AC 7 ^ = Jj_/3. 

The linear stability analysis of these mean-field solutions 
can be performed by looking at H16|l with p = 2 



aL(C 2 ) 



J±k 2 2 



2C 



1 



-7W 



(23) 



These eigenvalues must be positive, for all momentum k, 
in order to ensure the stability of the solution. 
Expanding first around the vertical line, C2 = 0, lll'-il) 
leads to the following stability criterion 

7 (k) = 7 log(l/fca) < 2, Vfc 



where 7d = 7(/ i3 1 ) and Id is the Debye screening length 
taken into account as a hard cut-off to eliminate the log- 
arithmic divergence of the Coulomb term in the thermo- 
dynamic limit. Details concerning the screening mech- 
anisms and the expressions of Id have been reported in 
Appendix A. As can be seen in the first term of 119|) 
in the confinement regime, the absolute value has been 
dropped. This is the non-return approximation for the 
confinement regime the limitations of which will be con- 
sidered at the end of this section. In the following sub- 
sections we will see how the tilted phase arises via zero- 
temperature phase transitions in both p = 1, 2 regimes. 



k > k c = - exp(-2/7o) 
a 

where k c is the Coulomb-dependent critical momen- 
tum. This shows that the instability is driven by low- 
momentum, i.e. long wavelength, modes. Reminding 
that 1/Id < k < 1/a, we find, in agreement with the 
mean- field analysis, that the vertical string is stable as 
long as k c < I /Id- k c = l/l D returns us to fflty. Increas- 
ing the strength of the Coulomb interaction increases k c . 
For k c > I/Id a continuous set of modes, i.e. between 
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FIG. 1: 1) Mean-field free energy of the elastic string 1191 as a 
function of £ exhibiting a second order phase transition. The 
two curves correspond to 7d < y D where the vertical string 
is stable and jd > Jr> where the tilted line is stable as can 
be seen from the double well. 2) Free energy of the confined 
string as a function of £ exhibiting a first order transition 
from the vertical string to the free string. From upper to 
lower curve: free model 713 = 0, metastable states appear 
when 715 > Jo, they become stable states when 713 > 7f>. 



1/Id and k c , has negative eigenvalues so that the vertical 
string is no more stable. 

Expanding around £2 7^ 0, we see from <|23|) . that all 
modes are stable as soon as £2 > 1/V2- For such tilt 
angles we are of course above the critical point as can be 
seen with the help of The tilted line is thus stable 

against harmonic fluctuations. 



hand configurations with £1 < l/y2 = £* are always 
unstable. This implies that the stability of the verti- 
cal string cannot be analyzed within the present model. 
This is related to the non-return approximation made in 
iJTSJ. Numerical simulations however show that the tilted 
string is, here also, the new ground state, cf. section 4. In 
relation with the crossover between elastic and confine- 
ment regime we can then interpret £* as a low boundary 
for the crossover angle Co, cf- ©■ I n the following we 
will thus consider solutions corresponding to £1 > £*, i- e - 

ID > 1* D - 

Figure n shows a crucial difference between the confine- 
ment and elastic regimes concerning the mechanism by 
which tilted strings, i.e. a double well in the free energy, 
appear. Contrary to the previous case the second deriva- 
tive of the energy remains positive a signature of the 
fact that the correlation length remains finite and that 
the transition is first order in the confinement regime. 
Metastable states thus appear. We are going to show 
first that they appear at 7^, the upper spinodal strength. 
From l|19|l with p = 1, the vanishing of the first derivative 
of the free energy density leads to 



1=7/3 



(1 + Ci 2 ) 3/2 ' 



(25) 



The spinodal line, above which (metastable) solutions 
appear can be defined with the help of the coupled equa- 
tions and it's derivative. This leads to 



1 

71 



Id 



3V3 



2.6 



(26) 



Increasing the Coulomb interaction we reach the criti- 
cal regime where the mean-field free energy of the tilted 
solution becomes equal to that of the £ = solution 
/ x (£ = 0) = /i(Ci). This equation together with 
define the critical point and lead to 




Yd « 3.3 (27) 



As in critical phenomena, for 7^5 > J D , the system jumps 
to the non-zero tilt angle solution corresponding to the 
absolute minimum of the energy. 



C. The charged string in the confinement regime 

When £ becomes larger than the crossover value £ > £0, 
cf. JSJ, the elastic approximation is no more valid and 
we reach the peculiar confinement regime with p = 1. In 
this case the eigenvalues (|16[) are given by 

A 2 , fe (£i) = J± 2Cl2 ~ 1 5 fc 2 7 (fc). (24) 

[i + C 2 ] 5 

We see from l|24[) that configurations with £1 > l/y/2 = 
£f are always stable against harmonic fluctuations. 
Moreover this is a pure Coulomb stability. On the other 



D. Quantum fluctuations 

At T = the string is submitted to quantum fluctu- 
ations. The latter give a contribution to the difference 
correlation function of the form 

< 8x(y) 2 >« ft logy, 

where y is along the string. Assuming that the distance 
between lines is of the order of Id, collisions between 
lines will take place on a scale y w exp^u due to these 
transverse fluctuations. These collisions are thus present 
on scales much larger than our upper cut-off Id and will 
therefore not affect previous results. 
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IV. THE CHARGED STRING AT FINITE 
TEMPERATURE 

We seek, in the present section, for the effect of temper- 
ature. In the following we consider the low temperature 
regime where the string needs to be quantized. From the 
expression of the full free energy, thermodynamic quan- 
tities can be evaluated. We give the expression of the 
heat capacitance. In the presence of thermal fluctua- 
tions, the angular degeneracy gives rise to defects con- 
necting the ground states of the string: the angular kinks. 
The N— kink problem is defined and the energies of the 
1— and 2— kink string configurations are given. Single 
angular kinks are activated. Bi-angular kinks are non- 
topological and are subject to a Coulomb-confinement; 
these excitations correspond to a Coulomb-roughening 
of the string in addition to the usual roughening. 



A. Thermodynamic quantities 

Apart from a logarithmic correction the spectrum l|16|) is 
that of phonons, a feature of the saddle point approxima- 
tion. We thus have a harmonic oscillator problem which 
is easily quantized 

E {nk} =£ (0) +^2(n k + hu ktP 



which is linear in T with an angle dependent coefficient. 
Such an expression can be usefully compared with exper- 
imental work on related physical systems. 




a) 1 \ b) 





where the units has been chosen such that uik,p has the 
dimension of an energy. 
This yields 



FIG. 2: Angular kinks shapes, a) Exact kink shape, b) Point 
like kink, c) Shape of the string with iV— kink excitations, d) 
Numerical result. 



fp = /f + T I ; g log (1 - exp(-^ p /T)) (28) 



where the second term is an entropic contribution and 
the zero point energy has been included in the first term. 
The general expression of the free energy density is then 



fp ~ ft 



(0) 



T 2 



12^.7]. [P(P-1)+ 7(01 



(29) 



where 



7(0 = 7 



2Q-1 

(f + C 2 ) 5/2 



and 70 is given by l(T%|l . 

It is simple, from (|29(1 . to compute the heat capacitance, 
at constant length, of the system. We obtain 



T 



6W Jj_ 



P(P - 1) + 7 



2C 2 -1 

(c 2 + 1) 5/2 



B. Excited states of the string 



The angular kink solution 



Coarse graining the system we consider the kink as a 
point-like excitation as in Figure ^jp. This amounts to 
neglect it's core on the scale of which elastic or confine- 
ment energy compete with the Coulomb energy to con- 
nect smoothly the two ground states as shown in Figure 
12^.. At finite temperatures a general configuration of the 
string is thus similar to Figure^ where an array of angu- 
lar kinks is present. This shape agrees with the numerical 
results displayed on Figure 01. 



(30) The energy of an array of A" kinks is given by 



7 



AE 



N 



N 

Wy_ 



y t,j=o J y 
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Vj+i 
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1 



^/(y-y') 2 + (x l+1 (y) 
I 



c j+i 



(y'W 



E 



(0) 
JV 



(31) 



where w y has been defined in (|18fl . J/jv+i = — 2/0 = +oo, 
the other yj's denoting the position of the defects, Xi(y) 
is a portion of string connecting point to point j/i 

and -E$^ is the Coulomb energy of the same string with- 
out kinks, i.e. at T = 0. The point-like kink approx- 
imation amounts to use the tilted ansatz ground state 
x%(y) = (-l) l C(2/ - Vi) + Xi(yi). Eq. (jnU is available for 
both elastic and confinement models. This is consistent 
with the fact that these models have similar long distance 
properties as shown in Appendix B with the help of a 
transfer matrix approach. However, even though exact 
results can be obtained with the help of (|31|) within the 
limits of the point-like kink approximation, the obtained 
expressions, in the general case, are quite intractable. We 
will therefore consider simple cases. 
First the 1— kink energy. The latter reads 



AE\ 



,v / l + C 5 



In 



1 + Vi + c 2 



(32) 



where L is the length of the string. The energy of the 
kink is lower than the total energy of the string which is 
over-extensive in the absence of screening, i.e. ~ LlogL. 
Nevertheless, Eq. (1321 shows that the kink has infinite 
energy in the thermodynamic limit, a special feature of 
the long range Coulomb interaction. In the presence of 
screening, the energy of the string becomes extensive, i.e. 
~ Llog/c and the single angular kink is activated 



iDWy 



■An 



a-yV 1 + C 

Thus, we have an exponentially small density of kinks 

AEi 



exp 



T 



(33) 



2. Transverse fluctuations 

Given the density of kinks (|33(l , the length of the string 
between two angular kinks is given by 



I 



1 

n 



Within this length thermal fluctuations also play an im- 
portant role. The difference correlation function reads, 
with the help of (|15|) , in the limit of large y 



< (5x(y) - 5x(0)) 2 > p « y 



Ta,, 



1 



Jl p(p-l) + 7(C)la(v/<0 
(34) 

where we have omitted higher order terms. Equation l|34(l 
shows that the charged string, which might be tilted if 
the Coulomb interaction is sufficiently strong, roughens 
at all non-zero temperatures. This roughening has two 
origins: there is a usual roughening due to the elasticity 
of the string and a Coulomb roughening, originating from 
the spontaneous symmetry breaking. Hence, there is no 
positional order 



< 5x >r 



y/< (5x(y) - 5x(0)f 



oo 



but there is long-range orientational order 



< 6£ >r 



d < 5x> 
dy 



Notice that the screened Coulomb interaction does not 
affect the value of the roughening exponent which is 1/2. 
This is also known from a theory of fluctuating charged 
manifolds, cf^, as has been said in the Introduction, and 
confirms the validity of our approach. What is peculiar to 
our system is the Coulomb roughening of the string due 
to the spontaneous symmetry breaking. The latter might 
be isolated from the usual roughening by the following 
procedure. 

We consider the 2— kink problem on the same line as 
the single kink of the previous paragraph. When their 
separation Sy, along the axis of preferable orientation, is 
smaller than the screening length, equation 131fl yields 



AE 2 = 



y/(l + C 2 ) P + 5y) 2 + C 2 (L - eSy) 2 } + L + dy + ( 2 (L - eSy) 



2Sy [1 + C 



8 



+Sy log 



y/Q. + £ 2 ) [L 2 + ( 2 (L-Sy- eSyT] + L + ( 2 (L-5y- eSy) 
26y [l + C 2 ^] 



+L\og 



+ C 2 ) P + Sy) 2 + ( 2 (L~ gffij + L + Sy + ( 2 (Sy - eL) 
2L[l + ( 2 ^] 



(L - (fy)log 



y/(l + C 2 ) [L 2 + £ 2 (£ -Sy- eSyf] +L + C 2 (6y- e(L - Sy)) 
2(L-Sy) [l + C 21 ^] 



(35) 



The term with e = +1 corresponds to the contribution 
of the kinky string whereas the one with e = — I is the 
subtraction of the ground-state string. When Sy < Id <C 
L, iJSSJ yields 

A^ 2 «^=l0g(l + C 2 ) + 0((^) 2 ). (36) 

In the limit L — > oo the angular bikink has a finite en- 
ergy. At the optimal angle, equation l|36|) shows that 
this pair is confined with an energy growing linearly as 
a function of the transversal distance between the kinks. 
The proliferation of angular bikinks is therefore shown to 
be a signature of the Coulomb roughening of the charged 
string. 

V. NUMERICAL APPROACH 

The statistical properties of the charged string in both 
elastic and confinement regimes has been studied with 
the help of a Monte Carlo-Metropolis algorithm. These 
numerical simulations allow the computation of the prop- 
erties of the full Hamiltonian IjlUI) with open boundary 
conditions as above. The ground state is found with the 
help of a simulated annealing. The shape of the distribu- 
tion, the value of the zero temperature order parameter 
£ and the energy of the distribution as a function of 7 
have been computed to detect the various instabilities 
and associated phase transitions. The effect of thermal 
fluctuations, i.e. the roughening and tunneling, have also 
been studied. These numerical results show good agree- 
ment with the analytical results presented above. 

A. Elastic string 

We focus first on the zero-temperature properties of 
the string. The ground-state has been found by per- 
forming a logarithmic simulated annealing during each 
Monte-Carlo simulation 

r n In 2 



where i is the Monte-Carlo step and T[0] is a high ini- 
tial temperature. Due to an improved algorithm to treat 
the Coulomb interaction the logarithmic decay of tem- 
perature has been successfully implemented in the sim- 
ulations. As is well known this leads to the ground- 
state with the lowest probability of being trapped in 
metastable states. 

We compute the properties of the full Hamiltonian i|10[l . 
The numerical coupling constant in @ is defined by 
Inum- The rigidity J± in J§J is taken to be unity. 
With these numerical coupling constants a comparison 
can be made with the analytical results obtained using 
a mean-field theory. In this respect 7„ urn should cor- 
respond to 7 in (|17fl . The computational cell needs of 
course to be finite and there is therefore no reason for 
using an upper cut-off limit Ijj. Therefore Id — » L and 
Inum = InumD / hi L which reflects the non-extensivity 
arising from the long-range interaction. In particular the 
critical numerical coupling constant should be given by 
Inum = InumD / lnL where Inum should be obtained by 
simulations and 7£ um £> compared with (|20|l . 

We consider first the behavior of the order parameter £ as 
a function of j num for different sizes L. A typical result 
is displayed on figure |3J Each point in this figure cor- 
responds to an independent Monte-Carlo simulation, i.e. 
given by generating a new random configuration which is 
then thermalized by the simulated annealing and brought 
to an effective zero temperature. Due to the double de- 
generacy of the ground state it is the absolute value of £ 
which has been plotted. The behavior is the one expected 
by the analytical results. Further more the critical cou- 
pling constant given by 7|j num agrees with 7^ with an 
error of the order of 1% as is given by fitting the numeri- 
cal data with (|2f I) . Below the transition the mean angle is 
zero. A typical shape of the string in this region is given 
in Figure The string is vertical in average. Above the 
transition the mean angle increases with j num in accor- 
dance with JUJ. The string is then tilted with respect to 
the main axis as shown also in Figure^] No jump in the 
order parameter takes place which is an indication that 
the phase transition is second order. The energy of the 



9 



Elastic regime, L=100 
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FIG. 3: Absolute value of £ as a function of 7, 
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FIG. 4: Shape of the string: below the transition the charged 
string is vertical. Above the transition it tilts away from the 
anisotropy axis. 

string has however been computed as a function of 7„ um . 
It is continuous. The second derivative is found on the 
other hand to diverge around the critical value. 

Next we take into account thermal fluctuations. When 
the temperature is less than the activation energy of sin- 
gle angular kinks AEi , cf . \.\'2\) , the line is rough as shown 
by Figure[S| Above AEi single angular kinks are present 
in addition to the roughening, cf. Figure [2 The orien- 
tational order is then lost. These results clearly confirm 
what has been said in the previous sections. At tempera- 
tures higher than AEi , thermal fluctuations average the 
tilt angle to zero. 



B. Confinement regime 

The same approach is performed in the confinement 
regime. 

Figure displays the absolute value of £ as a function 
of 7„ um . The behavior is the same as in the previous 




FIG. 5: Rough inclined string at low temperature; orienta- 
tional order is preserved. 



confinement regime, L=100 
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FIG. 6: Absolute value of ( as a function of 7„ um in the 
confinement regime. 



case. The straight string is stable below the transition 
and tilts above. However a jump in the order parameter 
is noticeable (especially in comparison with the equiv- 
alent Figure |3| concerning the elastic regime) , an indi- 
cation that the transition is first order. This is further 
confirmed by looking at the energy displayed on Figure 
It is clearly discontinues as a function of 7, mm at the 
transition. The numerical study allows to capture the 
properties of the charged confined string in the whole 
range of 7„« m . This is to be contrasted with the analyti- 
cal work where we could not consider the small tilt angle 
region. In contrast with what has been said inii this 
regime is not manifested by an instability of the string 
towards a modulation. The ground state given by the 
numerical simulations is also the straight, either inclined 
or not, string. 

At finite temperature the results for the confinement 
regime are qualitatively the same as for the elastic one. 
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confinement regime, L=1 00 




FIG. 7: Energy of the string as a function of 'Jnum- The solid 
line corresponds to the energy of a string. A discontinuity 
is clearly noticeable. The diamonds correspond to the first 
derivative of the latter. 



VI. APPLICATIONS 



A. The solitonic lattice (2D) 

The present study should be connected, as has been 
said in the Introduction, to the previous work on the 
statistical properties of solitons. Above, we have been 
considering the domain-line as given. Actually, in quasi- 
one dimensional systems, it has been shown inii, with 
the help of a mapping to the ferromagnetic Ising model, 
that there are no strings, i.e. infinite domain lines of 
solitons, at non-zero temperature. The latter emerge 
from a process of exponential growth of aggregates of 
bisolitons below a crossover temperature T cro . In the 
non-interacting case, T cro = Tq and the aggregates are 
finite rods perpendicular to the chains of the quasi-one 
dimensional system. At zero temperature the rods cross 
the whole system forming the free domain lines of soli- 
tons. In the interacting case the crossover temperature 
is lowered with respect to T and becomes Coulomb de- 
pendent. In the region between To and T cro , intrachain 
antifcrromagnetic interactions, within the ferromagnetic 
ground-state, take place in order for the bisolitons to ag- 
gregate, cfiii. This is due to the important size-charging 
energy of the aggregates. A crucial approximation is that 
we have been neglecting shape instabilities in the inter- 
acting case. This reduced us to the region of the phase 
diagram where jd < 7£>, c f- P0|> - At high coupling 
constants, jd ^> j D , strings at T = 0, or aggregates at 
finite temperatures should melt into a Wigner crystal of 
solitons. At intermediate coupling constants the present 
theory, dealing with shape instabilities should be rele- 
vant. It should be noted that in the previous paper we 
have suggested that shape instabilities give rise to a mod- 
ulated string which would correspond to a single negative 



mode; as we have shown, the instability at intermediate 
coupling constant, is towards a tilt of the interface and 
affects a continuous set of modes. 

We will first seek the stability of the string phase with 
respect to the Wigner crystal of solitons at T — 0. Then 
we will concentrate on the non-zero temperature regime 
and look for the shape of the aggregates of solitons which 
lead to strings. To make comparison with the previous 
paper easier, notice that the present coupling constant, 
7 — Id I log Id , is similar to the coupling constant used 
in the previous paper T. Otherwise we have tried to 
preserve the same notations. 

1. Zero temperature results: from strings to Wigner crystal 
of solitons 

We define 7^j c as the coupling constant above which 
strings melt into a Wigner crystal at T = 0. We would 
like to check here that 7$° > 7£> in order to preserve the 
coherence of the arguments. 

As the Coulomb interaction increases, the string of soli- 
tons tilts. It's energy also increases. This leads to a 
dissociation of the string as soon as it's energy per soli- 
ton [i = Silt ' fSN becomes larger than the energy of 
the constituent soliton n so i. The latter corresponds to 
2J±R/a x — (ze) 2 /eR where R » \J sjv is the radius of 
the Wigner Seitz cell including the soliton and a back- 
ground charge of density v/s with s — a x a y . This energy 
can also be written as 

2J X R h , 
V>sol = [1 - V71 

where 7 is given by iJTSJ. The cohesion energy of the 
string is thus 



E c = Jx\C + 



Id 



(37) 



From Q37[l we see clearly that in the absence of Coulomb 
interaction, which implies that £ = 0, the vertical string 
is stable. When the Coulomb interaction increases the 
tilt angle increases beyond the transition to the tilted 
phase. Thus E c increases from negative values. The 
melting of the string is given by E c = and, Q being the 
optimal angle of the string, l|25|l. Starting with 7d«1 
in E c — 0, leads to the boundary 723 = l/v ^> 1, 
in contradiction with the starting hypothesis. On the 
other hand, starting with 7^ ss \ jv leads to the lowest 
boundary 70 w j Dl cf. (|26[) . again in contradiction with 
the starting hypothesis. The solution satisfying the self- 
consistent equation is thus in the intermediate coupling 
constant range j D < 7^ c <\/v and is given by 

Id « 1/4^- (38) 
We see from (I38|l . that for a sufficiently diluted system, 
7d c > 7d> cf - Tne tilted phase might thus be 

observed when this condition is satisfied. 
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FIG. 8: a) Rod in the absence of Coulomb interaction. The 
Coulomb energy of the lozenge b) is lower than that of c) . At 
finite temperature the aggregates thus have the b) shape. 



2. Non-zero temperature results: shape of the aggregates 

More exotic is the non zero temperature case. From the 
point of view of interfaces, the latter roughen and thus 
collide at non-zero-temperatures. Equivalently, from the 
point of view of solitons there are no infinite lines but fi- 
nite size aggregates. The size of the aggregates is approxi- 
mately given by the distance between two successive colli- 
sion points of neighboring interfaces. We address here the 
question of the shape of the charged aggregates of soli- 
tons in the intermediate regime where 7f, < 7d <C Jr> c - 
As has been said above, in this regime, the aggregates 
must adapt their shapes in order to minimize their elec- 
trostatic energy. An inclined aggregate, as in Figure |H| 
c), would thus be more favorable than the straight rod. 
But at finite temperature angular kinks are present and 
they further reduce the Coulomb energy. We can shown 
that a droplet with the lozenge shape of Figure[S]b) is fa- 
vored in comparison with the previous aggregates. This 
is done by evaluating the energies of these aggregates 
AE b = E b - E ( h se) and AE C = E c - E ( c se) where their 



divergent self-energies E^^ and E^' have been sub- 
tracted. For the aggregate c) the separation between the 
bikinks x op t is taken as the equilibrium separation under 
the combined action of the Coulomb force (ze) 2 /ex 2 and 
the confinement one —J±/a x , i.e. x op t = y/sj- Moreover, 
the tilt angle corresponds to the optimal one. In the con- 
finement regime 7 is related to ( with the help of ()25ll . 
The regularized energies of the aggregates thus depend on 
one free parameter (. The later must be larger than 1 / \[2 
to be in the stable confinement regime. We find that the 
difference between these energies AE C — AE b > as a 
function of £ so that the energy of b) is always lower than 
that of c) above the critical coupling constant 7^, . 



? (se) 



The results obtained in this section are summarized on 
the phase diagram Figure 
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FIG. 9: Phase diagram of the 2D system as a function of the 
coupling constant 70 and temperature T. Details are given in 
the text. 



B. The solitonic lattice (3D) 

We consider briefly the important case of three dimen- 
sional systems where interfaces are domain walls. In 3D 
some rigorous statements allow a good understanding of 
the uncharged case in the frame of quasi-one dimensional 
systems, cf4&. In the presence of Coulomb interactions 
we shall thus follow^ as well as our present experience of 
the charged 2D case. Concerning notations, in the third 
direction, 2, the coupling constant will be taken as J z 
and the unit length as a z . 



1. Uncharged domain wall 

For the non-interacting case it has been shown inifl that 
peculiarities are brought by the raise of dimensionality. 
In 3D the excitations analogous to rods would be disc-like 
objects. However, rigorous statements concerning the 3D 
case show that, at a temperature T3 « 8J±/ ln(2/i^) <C 
T c « Ji/v, a density v c = v — v n of solitons condense 
into infinite anti-phase domain walls perpendicular to the 
chains of the quasi-one dimensional system; v n then cor- 
responds to the density of finite size solitons below T 3 . 
Therefore, the crossover regime of growing rods, which 
took place below To in 2D, cannot take place in 3D and 
is replaced by the transition at T3. 

The effect of thermal fluctuations on the walls is well 
known. Several authors have shown (see references in™) 
that, in 3D the hight-hight correlation function remains 
finite at T < Tr whereas it diverges logarithmically for 
T > Tr, where Tr > is the roughening temperature. 
From these considerations T3 may be considered as the 
roughening temperature. The existence of walls in 3D 
at finite temperatures is compatible with the fact that 
in Ising-like systems there is 2D long-range order at low 
temperatures. 
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2. Charged domain wall versus Wigner crystal 

The three-dimensional case of charged walls can be 
treaded the same way as the two-dimensional case. We 
will therefore consider again the elastic and confinement 
regimes. The Coulomb interaction will lead to a tilt- 



J 



ing of the wall here too. Two tilt angles can be defined 
now, because the anisotropy axis becomes an anisotropy 
plane. We will thus consider ( y which corresponds to the 
previous, in plane, tilt angle and Cz which is the tilt an- 
gle with respect to the third direction. The mean-field 
hamiltonian of the system is then 



H P /Jy = C? + axC? 3 



Id 



- 1 



(39) 



where, a± = J z /J y , 7d = e 2 ^ / 'es±J y , s± — a y a z and in 
the last term, the contribution of the straight domain wall 
has been substracted. Notice that it is only in the case 
where Cy — Cz = C that we gain a global factor 1 / y/l + ( 2 
as in the 2D case. The logarithm is therefore the signa- 
ture of a non-homogeneous charge dilatation within the 
wall. 

Our aim is to determine the phase diagram, jd versus a± 
in both regimes, i.e. p y>z — 1,2. Among the four possible 
cases we will consider only the full elastic (p y — p z = 2) 
and mixed (p y = 2, p z = 1) regimes. 

In the elastic regime, p y — p z — 2, the phase diagram 
is displayed on Figure ITU1 For 725 > 2a±, the stable 
solution is given by 



Cy = Cz = 



(40) 



In the region 7^ < 2a±, the stable solution corresponds 
to 



Cy = C* 



Id 
2a 1 



(41) 



For Yd « 16 < Id < (2oa) 3/2 ^/logoa /4 a third phase, 
corresponding to 



Cy 



(?'-(?)) 



1/3 



Cz = o 



(42) 



is stable. The domain of validity of this solution overlaps 
the ones of the two previous solutions. With the help 
of energetic considerations, we find that a line of first 
order critical points, at 7^ 20, separates phase l|40f) 
from phase (|42l) . By the same way, a line of first order 
critical points has been found to separate phase (|12fl from 
phase (|42|) . The line jd = 2a ± is a line of second order 
critical points which terminates at the critical end-point 
a± = 10, 7d = 20. Phases which would have both angles 
none-zero appeared to be unstable. 

A similar analysis can be done for quasi-two dimensional 
(Q2D) systems where the interface is elastic within the 
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FIG. 10: Phase diagram of the elastic wall: 70 as a function 
of a±. 



planes, i.e. p y = 2 and in the confinement regime between 
planes, i.e. p z = 1. For Q2D systems, the coupling in 
the third direction is weak with respect to the coupling 
in other directions, i.e. a±_ < 1. In the region of the 
phase diagram where a±_ < 1 the solution is thus given 

by 



Cy 



Cz 




(43) 



and is stable for 70 > 2a±. As before, when jd < 2a±, 
the untiltcd solution is stable for small a±_ and is sepa- 
rated from the previous by a line of second order critical 
points. From (|41|l and Ij43(l . we see that for a given value 
of the ratio 7d/2q;_l, the tilt is stronger in the mixed 
regime than in the pure elastic regime. This can be inter- 
preted by considering the equilibrium position of a con- 
stituent of the wall due to the constant confinement force 
and the constant electric force due to the charged wall; 
a priori, the equilibrium is reached only at 7^ = a 
feature which is peculiar to the confinement regime. This 
implies that the wall should melt as soon as > cx± or, 
if the Coulomb energy scale is sufficiently weak, that the 



13 



wall should strongly incline to reduce the coupling con- 
stant. Allowing the wall to tilt, the equilibrium equation, 
7d = Q!j_, is generalized to 




which leads to the following minimal stability angle 

t - .,(,(. -£))-!. 
above which the wall stabilizes at the equilibrium angle 

In 3D, the situation is thus richer than in the 2D case 
with the appearance of a tricritical point and quantita- 
tive differences between the various regimes but, qualita- 
tively, the same phase diagram for all of them. In partic- 
ular, the tilting of the interface takes place with respect 
to a single anisotropy axis and not both. In quasi-two 
dimensional systems, the wall should tilt with respect 
to the third direction rather than y— direction and the 
tilting is stronger in the confinement regime to prevent 
the melting. Another difference with the 2D case is that 
thermal fluctuations do not affect these results, as long 
as we are below the roughening temperature of the wall. 



C. Stripes in oxides 

Up to know the whole theory is model independent. As 
we have already mentioned our results could be applied 
to doubly commensurate charge-density waves such as 
polyacetylene. We will however consider physical systems 
which are quite extensively studied nowadays. These are 
the oxides, e.g. cuprate, nickelate, manganese and the 
new ferroelectric state. The latter will be considered in 
the next section. We will focus here on the phenomenol- 
ogy of stripes with the help of the above theory. We con- 
centrate only on the long distance properties, namely, the 
fact that stripes consist of lines or walls of holons which 
are similar to the charged solitons we have been consid- 
ering. The latter are thus charged interfaces as the ones 
described in this paper. 

The oxides we consider are quasi-two dimensional materi- 
als formed by coupled planes. At intermediate values of 
doping, a 31? ordering of the stripes might have been 
seen in La\^- X N do.4Sr x CuOi and La2- x Sr x Ni0^j r s, 
cf£ and references therein. The impact of the pinning 
by the lattice and the Coulomb interactions play a cru- 
cial role in these systems. 

In Lai,Q^ x Ndo.4Sr x Cu04 neutron diffraction studies 
have revealed diffuse magnetic peaks along the z- 
direction, in the stripe phase, which is consistent with 
stripes as one-dimensional interfaces. Moreover, X-ray 



studies reveal that stripes are parallel in next-nearest- 
neighboring layers but orthogonal between nearest lay- 
ers. This 3D ordering is probably due to a pinning of the 
stripes by the lattice structure. 

The case of nickelates is more interesting. In 
La2- x Sr x NiOi+s there is a body-centered stacking of 
parallel elastic stripes with only weak perturbations due 
to the lattice. This 3D ordering is due to Coulomb re- 
pulsion and might correspond to a domain wall with a 
step-like profile suggesting a confinement regime in the 
third direction. This justifies the mixed confinement- 
elastic regime we have considered in IVIBI Experimental 
results do not however exclude the interpretation of this 
3D ordering as a Wigner crystal of lines. The latter hy- 
pothesis allows us to focus on strings within planes of the 
quasi-two dimensional system, as we are going to do in 
the next paragraph. 

Following Wakimoto & all, we are going to recall briefly 
here some experimental results related to the tilting of 
string-like stripes in Lai- x Sr x CuO \. 
The experimental results are summarized on Figure 1111 
which has been taken froml. In Figure ITU (a). 5 = 2n/d 
where d is the distance between stripes. At low doping, 
the inset on the left of this figure, shows magnetic Bragg 
peaks corresponding to a one-dimensional spin modula- 
tion which is inclined with respect to the reference tetrag- 
onal axis. At higher doping, the inset on the left, shows 
that this spin modulation, along with the charge stripes, 
have rotated by an angle of 45 degrees. The main fig- 
ure, together with Figure ^TJ (b) which displays the tilt 
angle as a function of the concentration of dopant, show 
that this shape instability takes place around x c ~ 0.05 
where 8 increases linearly with x. Moreover, the authors 
report a weak dependence of the magnetic peaks on the 
third direction. This implies that the spin modulations 
are weakly correlated between successive C11O2 layers. 
Thus, the stripes can be considered with a good accuracy 
as charged strings which justifies our starting hypothesis. 
The stripe instability observed experimentally might be 
due to the competition between the long-range Coulomb 
interaction and an anisotropic interaction. To relate 
this instability to an increase of the concentration of 
dopants x we should probably take into account struc- 
tural changes. As argued by Tranquada & a£, the lines 
of charge might be pinned along the direction of the tilt 
of the CuOq octahedron. This tilt direction varies with 
doping. It is along the orthorhombic axis at low doping 
and along the tetragonal axis at higher doping. These 
axis are 45 degrees from each-other. At an intermedi- 
ate doping, corresponding to the critical value x c , a 45 
degrees tilt of the stripes, satisfying the structural con- 
straints, can take place. 

Simple considerations related to the competing energy 
scales are in favor of such an order of magnitude for the 
tilt angle. First, in this system e 2 /ea ~ 0.1 — XeV where 
a is a unit length of the tetragonal structure. Second the 
exchange energy J± is of the order of Q.13eV. Thus 7 ~ 
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FIG. 11: Tilting of stripes in a cuprate oxide (from—). 

0.4—4.5 where (|18f) has been used. If we consider that the 
Debye length is of the order of the inter-stripe distance 
1/(5, at the transition x c we have Id ~ 20a. Thus 7d = 
jlnQu/a) ~ 1.2 — 13.5. Using (|21[1 . the maximum tilt 
angle compatible with these values is around 55 degrees. 
A tilt angle of 45 degrees can be reached with a 7£> ~ 4.0 
which is clearly in the range of allowed values above. 
Hence, the competing interactions can be responsible for 
such a tilt. 

This tilt could reduce the electrostatic energy of the 
stripe as can be seen by the following argument. At low 
doping, i.e. in the diagonal stripe phase, the hole con- 
centration along the stripe is 0.7 hole/Cu. This concen- 
tration is reduced to 0.5 hole/Cu in the collinear stripe 
phase, cf. Wakimoto & al. The density of charges along 
the stripe is thus reduced by y/l + £ 2 , with £ = 1 for the 
present 45 degrees tilt angle. This dilatation factor re- 
duces the electrostatic energy of the string as has already 
been seen in <|19ll . This could hold assuming that the 
number of holes/stripe is constant. The added dopants 
would lead to the formation of new stripes. The latter 
become tightly packed at high doping as suggested by 
the linear relation 5 « x. 



D. Uniaxial ferroelectrics 

Recently a ferroelectric Mott-Hubbard state has been 
observed in quasi ID organic superconductors such as 
(TMTTF) 2 X . We give a brief account on the origin of 
this phase and connect the results to our present study. 
Experimental and theoretical results are given in detail 
irA 

In the system mentioned above TMTTF are molecules 
composing the stacks and therefore giving rise to the one 
dimensional nature of the compound. X are ions placed 
in the vicinity of the molecules. At large temperatures, 



of the order of lOOif, the ions reorder through a uni- 
form shift which gives rise to a ferroelectric phase. The 
latter has been detected via the measurement of the di- 
electric susceptibility which shows a gigantic anomaly at 
low frequency around Tf e ~ 150-fT. At low temperature 
domain lines, separating domains with opposite polariza- 
tion, should be observed. Their existence is necessary to 
minimize the external electric field generated by charges 
accumulating at the boundary of the sample. Moreover, 
as the interface corresponds to a jump in the polarization 
it is charged. Following^ we shall define as a— solitons 
the elementary constituents of the interface which con- 
nect domains with opposite polarization. They carry a 
fractional charge ze with z — 2a /ir. 
Two facts suggest that this system is a good candidate 
for the observation of charged interfaces. First the back- 
ground dielectric susceptibility eo ~ (uj p /A) 2 ~ 10 4 is 
large. Further more 7r— solitons are also present and 
they screen the Coulomb interaction. This weakens 
the Coulomb interaction. In such systems, the strings 
of a— solitons should thus be stable against a melting 
towards a Wigner crystal. Bubble-like aggregates of 
a— solitons, as the ones described in IVII could also be 
observed at higher temperatures. 



VII. CONCLUSION 

The statistical properties of charged interfaces, strings 
and walls, has been studied. We have found that shape 
instabilities, due to competing interactions, play a fun- 
damental role. They manifest themselves through zero- 
temperature phase transitions from which new ground- 
states emerge where the interface is tilted. The 2D case 
of the string has been extensively studied. For small tilt 
angles the string is elastic and the transition to the tilted 
ground state is second order. On the other hand a first 
order transition takes place if we are beyond the limit 
of validity of the clastic regime, i.e. in the confinement 
regime. At non-zero temperatures it has been shown 
that in either regimes a pure Coulomb roughening of the 
string, due to the proliferation of angular-bikinks, exists 
in addition to the usual roughening as soon as T > 0. 
At higher temperatures angular-kinks are thermally ac- 
tivated and connect the degenerate ground-states. These 
results have been confirmed with the help of a numerical 
approach based on the Monte Carlo - Metropolis algo- 
rithm. We have then related the present study to the 
general problem of the statistical properties of solitons in 
charge-density wave systems in either 2 or 3 dimensions. 
In 2D systems the string has been shown to emerge from 
lozenge-like aggregates of solitons upon decreasing tem- 
perature. In 3D, the additional coupling constant in the 
third direction enriches the tilted phases with the appear- 
ance of a tricritical point separating them. In the confine- 
ment regime, a disintegration of the wall is favored but 
the solitons might order in the form of a strongly inclined 
wall. Applications concerning stripes in oxides and a pos- 
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sible explanation of their tilt, as well as charged strings 
in uniaxial ferroelectric experiments, have also been con- 
sidered. 
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Appendix A 

In this appendix we give some details on the screening in 
the ensemble of solitons. Screening makes the energy ex- 
tensive and gives a meaning to the thermodynamic limit 
as has been said all along the paper. In the following we 
will work in the RPA approximation which is very well 
known in the field on many-body theory. This allows us 
to derive the Debye length and the electrostatic potential 
in two cases: two and three-dimensional screening. 
We consider first a single plane embedded in 3D space. 
This case corresponds to two-dimensional screening. 
Supposing the Coulomb field <f> slowly varying, it obeys 
the following semi-classical equation 

-eA(j)(R) = -47re 2 (n (0) * <j)){R)5(z) + A-ne 2 P Q{~r)8{z) 

(44) 

where R is a 3D vector associated with the 3D Coulomb 
field, V is a 2D vector indexing the charges which are 
embedded in the z = plane and po is some external 
charge. IT ) (x, y) is the response function, or polariza- 
tion part, of the soliton system which is related to the 
correlation function by fluctuation dissipation theorem 



dispersion via the convolution. The Coulomb field is then 



TH {0 \~r) =< 6p(r)6p{0)> 



(45) 



where 5p(r ) — p(r) — p(r) being the density of 
solitons and v the background neutralizing charge. 
The first term on the right of (|44|) corresponds to the 
screening charge. This term implies linear screening and 



[k ,z = 0) = 



2we 2 po(k) 
:\k\ +27re 2 lK )(~fc ,z = 0) 



(46) 



where k is the reciprocal vector associated with r . The 
fact that the potential is evaluated at z = implies that 
the field lines are confined to the same space as the parti- 
cles, cfpii. Expression (|46|l leads to the Debye screening 
of the field with a Debye length 



1 



2ne 



— = _n<°>(fc = o). 



(47) 



Also the k— term in the denominator of l|46|) is a feature 
of 2D— systems. In real space this leads to an algebraic 
screening instead of the exponential one, cf^£ for detailed 
studies of screening effects in Coulomb fluids and below 
for an example. 

We consider now the three dimensional case of indepen- 
dent planes embedded in a neutral media. Such a system 
of decoupled planes brings a 3D— screening. The screen- 
ing charge is then a function of the three-dimensional vec- 
tor R. But, as the planes are independent, the bare cor- 
relation function corresponds to H^(R) = H(°\~ : r)S(z). 
By the same arguments as in the first section we can com- 
pute the coulomb field which reads 



4ire 2 p Q (K) 



eK 2 + 47re 2 n(°)( k 



0)/o, 



(48) 



where we consider only the z = plane subject to the 
bulk screening. The three-dimensional Debye length is 
then given by 



1 



Ane 2 



I 2 

L 3D 



n(°)(fc = o) 



(49) 



where a z is the inter-plane distance. In real space l|48|l 
leads to the usual exponential screening. 
Finally we derive the polarization part for the present 
problem. This gives explicit expressions for the screening 
lengths and for (p(r ). Below the crossover transition To 



TH^O, y) w ^ua x exp (-^- j (2S(x) + S(x + a x ) + S(x - a x )) 



(50) 



ul. 



corresponding to the aggregation of the solitons into rods Thus D^ ) ( k ) ~ ., , 
of length l x , dM. Fourier transforming |JHJ we get gity of solitons v 



and, as we work at constant den- 
\/2^ exp(2 J±/T). With the help 



4sa 



v , , I 

1 + cos{k x a x )) 



16 



of |g7J| and ijlHl, this yields 



2esa y T 



TIC 



! (2i/) 3 /2 



exp 



'3i) = 



/ esa z T 
7re 2 (2i/) 3 / 2 



exp 



T 



(51) 



In the 2D— screening regime, the field can then be writ- 
ten, in the long distance limit with respect to the coarse 
graining length a x 



2ire 2 



1 



k v l ± 



/4 



fc 1 



(52) 



which reflects the anisotropy of the problem. To give an 
idea of the algebraic screening we Fourier transform 1)521) 
in the limit k y l± — > 0, i.e. far enough from the rod. This 
yields 

_ _ 27TC 2 1% 

in the limit r ^> Id- We see that the algebraic tail leads 
to a slow decay of the potential in contrast to the usual 
exponential one. However this type of interaction, of the 
dipolar-type, beyond the screening length are compatible 
with the growth of the aggregates, cfiii. For more details 
on the peculiar screening properties of such system with 
anisotropic intrinsic charges see 



U9 



Appendix B 

Our purpose is mainly to show the differences and com- 
mon points between the two interface models. The ap- 
proach is based on the transfer matrix technic. In the 
present case, this formulation cannot be used straight 



forwardly because of the long-range interactions in 11011 . 
Some results obtained in the non-interacting case should 
however remain in the presence of the long range inter- 
action. It is on such a property that we focus now. 
The associated Hamiltonian of the problem is 



K 



L-l 



y=o 



(53) 



where the short distance cut off has been taken explicitly 
into account. The transfer matrix T a is defined as an 
integral operator with eigenfunctions il>( x ) whose kernel 
is expressed in terms of (I53f) . Namely 



K 



\x-x'\ p VO')- (54) 



The free problem is simple enough to be treated just 
with the integral equation l|54)l . For a more exhaustive 
account on the general results one can obtain with help 
of the transfer matrix technic seei&. 
For both models the eigenfunctions of the transfer matrix 
are simply plane waves 

ipk( x ) ~ exp(ikx). 

This amounts to Fourier transform the kernel of 154(1 . 
The eigenvalues of the transfer matrix associated to both 
interface models are thus 



2K 



„sos 



log (K 2 + k 2 ) 



(55) 



By comparing the two spectrums in (|55|l we see clearly 
that both models have the same long distance, i.e. k — > 0, 
properties. 
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